NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



^ ; o ^ 

(NASA-CR-164335) SYNCHROTEON EMISSIVITY N81-23868 

1 FROH MILDLY RELATIVISTIC PARTICLES (Stanford 
Univ.) 35 p HC A03/MF A01 CSCL 20H 

I Unclas 

G3/72 42359 i , 


synchrotron emissivity from mildly 

RELATIVISTIC PARTICLES 


by 

Vahe Petrosian 


I 



1 

! 


SUIPR Report No. 832 
April 1981 


National Aeronautics and Space Administration 
Grant NSG 7092 


\ 


W 




SYNCHROTRON EMISSIVITY FROM MILDLY 
RELATIVISTIC PARTICLES 


by 

Vahe Petrosian 


SUIPR Report No. 832 


April 1981 


National Aeronautics and Space Administration 
Grant NSG 7092 


Institute for Plasma Research 
Stanford University 
Stanford, California 



SYNCHROTRON EMISSIVITY FROM MILDLY 
RELATIVISTIC PARTICLES 

Vahe Petrosian 

Department of Applied Physics 
and 

Institute for Plasma Research 
Stanford University 
Stanford, CA 

ABSTRACT 

We present approximate analytic expressions for evaluation of the 
frequency and angular dependence of synchrotron emissivlty from mildly 
relativistic particles with arbitrary energy spectrum and pitch angle distri- 
bution In a given magnetic field. Eqs. (11) to (14) give our result In Its 
most general form. We show that eq. (14) can be simplified considerably as 
In eq. (22) and that for most practical cases eq. (11) accepts simpler forms 
(eqs. 29 or 35). These results are valid for particle distribution which ar? 
not extremely anisotropic. 

Our result agrees with previous expressions for a non-relativistic Max- 
wellian particle distribution, and when extrapolated to non-relativistic and 
extreme relativistic regimes, It agrees with the previous expressions obtained 
under those limiting conditions. 

We compare the result from our analytic expression with results from 
detailed numerical evaluations. We find excellent agreement not only at fre- 
quencies large compared to the gyro- frequency but also at lower frequencies, 

1n fact, all the way down to the gyro-frequency, where the analytic approximations 
are expected to be less accurate. 


Subject headings : Synchrotron Radiation; Sun: Flares, Radio Radiation; 

Stars: Accretion, White Dwarfs; Plasmas 






2 


I. INTRODUCTION 

Synchrotron radiation from highly relativistic particles is an old hat worn 
by astrophysicists successfully for many decades. At extreme non-relativistic 
energies the primary contribution to emissivity is the well known cyclotron radia- 
tion at the first and perhaps few low harmonics. In the intermediate range 
there does not exist a simple formula for evaluaiton of the emissi\ itity (or 
absorption coefficient) from particles with an arbitrary energy spectrum and pitch 
angle distribution. The usual practice is to use lengthy numerical methods '^or 
evaluation of these quantities at mildly relativistic energies (cf .e.g., Samaty 1969 
and Takakura 1972, in connection with solar flare problems and Lamb and Master 1979, 
in connection with accreting white dwarfs). Earlier, in connection with fusion 
plasma, Trubnikov (1958) and Drummond and Rosenbluth (1963) derived equations 
for emissivity (or absorption coefficient) of synchrotron radiation from a 
thermal plasma (cf., e.g., Bekefi 1966, ch. 6, and references cited there). 

In most astrophysical problems the low harmonics of cyclotron radiation 
will be either self-absorbed, or absorbed by surrounding plasma, or suppressed 
by Razin-Tsytovich effect. Consequently, only radiation at high harmonics 
where optical depth ^ 1 is of interest. We find that at such high harmonics 
there exists a simple approximate method for evaluaiton of the spectrum and 
directivity of the synchrotron radiation from an ensemble of electrons in a 
given magnetic field. In this paper we present the results of this method for 
the total emissivity at high harmonics. In ill we present the general formula 
and the details of the derivation of the emissivity for an arbitrary energy 
spectrum and pitch angle distribution of particles. The emissivity for some 
commonly encountered particle distributions are presented in §III where we also 
give some emperical relation valid for a large range of the parameters and 
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compare the results from these analytic expressions with those obtained by 
detailed numerical calculations. A brief summary of the result 1s presented In 
§IV. Polarization and Intensity of radiation In optically thick regions and at 
lower harmonics will be discussed In a future paper. 
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II. ANALYSIS AND GENERAL RESULTS 

Consider particles with charge e, mass m and distribution f(v,Y) where 
fdudy Is the number density of particles In the energy' (1n units of me ) Interval 
Y to Y + dy and with pitch angle cosine between p to p + dp. The specific 
emissivlty j^(0) of the synchrotron radiation at frequency v and at an 
angle 0 with respect to the magnetic field of strength B is given by 
(cf. Bekefi 1966) 


where 


2ire^VL, /•« r<-l 

*— 7 — 2” / dY| dpf(Y,M)vn (e.Y.y). 


( 1 ) 


n^(e,Y,y) • |(Cose- 3 p)^J^(x) + (l- 0 pCos 0 )^[xj[„(x)/m]^ J 6(y). (2) 

y = mv|j/Y - v(l-3pcose), x = vYBsin 0 (l-ii 2 )^v^ , 

( 3 ) 

''b “ 2?mc ' 10®(B/gauss)Hz. 

2 

Here J (x) is the Bessel function of orfer m and = 1 - v" . 

As mentioned above, we are interested in emission at frequencies v such 
that optical depth s 1. Under most astrophysical circumstancer this occurs 
at frequencies much higher than the gyro frequency We, therefore, seek 

solution to eq. (1) for (v/vj^) » 1. In this case the delta function 
indicates that m is large and that the relevant harmonics are closely 
spaced so that the sunination can be replaced by an integration. This integral 
can be carried out easily with the help of the delta function. This step then 
eliminates the summation, introduces a multiplying factor of y/v^ and replaces 
m by yv(T-p3cos0)/v|j. 
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Furthermore, since m is large and we are not in the extreme relativistic 
regime, the Bessel function and its deviation can be approximated as follows 
(cf. e.g. , Abramowitz and Stegum 1970} 


0^(x) « (2iTm)'*^(1-z^)‘^s"’ , xJ|^,(x)/m « (1-z^) J^(mz), 


ze 


2 

(1-z^) 


2 

i+d-z^) 


_ X , Bsine 


(4) 


m 


l-uBcose 


3/2 


This approximation is valid if mil-z*") » 1. Note that z s 1 and that it 

_ 2-2 
has its largest value when 0 “ j ®f'd y » 1 .0, i ,e. , z < 3 or 1-z > y so that 

the above approximation is valid for y^ « Vvjj. At higher energies 
y » 1 and v/vjj » 1, one is in the ultrarelativistic regime where the 
spectral and angular distribution of the synchrotron radiation are well known 
(cf. i.g. Bekefi 1961, Ginzberg and Syrovatskii 1964). With the above approxi- 
mations eq. (1) becomes 


' C 


— ( f* V(9,y,h)Z.^”'(6.T.m) 


(5) 


Y(y.y.9) . I (1-z^ . )(1-|i . 6c . ose)^ ^ m - f (1-6ycose). 

(1-Z^) (1-M3COS0) ^ 


1) Integration Over Pitch Angle . In general, Y is a slowly varying 

Pm 

function of \i while, because of the large exponent, Z varies rapidly 


6 


from its maximum value at y * to zero at m « ±1. We, therefore, use 
the method of steepest descent to evaluate the integral over y. This gives 
m • (vY/v|j)(l-y^^$cose) and 


e Vi 

Jv<«) '■-r 


( 6 ) 
h 




where is the solution to the transcendental equation, 

max 


M - SCOs6(l-e,), c, • [inS- • (?) 


For slowly varying pitch angle dist;ributions dlnf(p,Y)/dy « v/vjj (see below 
for a more precise limit), the last term in the square brackets can be ignored. 
Then as we shall see below, is of the order of vj^/v « 1 so that 
* 0 COS 0 and eq. (6) becomes 

niaA 



{f-)^ f ” [0+2 COt^0/Y2)(l-B2cOS^8)*‘]2^^ 

■'l 


(8) 


where 


^max 


9 ^ 

t > 

1 + (i+t^)** 



(l+t^). 


t 2 0YSine 


(9) 


For 9 = J *(vin,ax “ ^ “ oj, the expression in the square bracket is 

equal to unity and 


®max (■'^/2 ,y, 


^max^ = 




( 10 ) 


in agreement with the result obtained by Trubnikov (1958). 
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2) Integration Over the Energy . In general, the expression in the square 

bracket of eq. (8) varies slowly with However, 2^"’ increases rapidly with 
particle kinetic energy [as (y - i)^"*] when y-1 « 1 and approaches unity 
as Y , while decreases rapidly (as a power law, f«(Y-l)"? 

or exponentially, f « with particle energy. Most of the contribution 

to the integral comes from the vicinity of the maximum of f2^”’/Y. Again, 
a good approximation to the integral is obtained by using the method of 
steepest descent which gives 

J^(8) • {.e\/c)(v/v^)\l+2 cot^e/Y2)(l-9W5)^f(Scos6,Y^)a^J^(t^)X. 01 ) 

where m * (VYgV^lCl+t^ ) and y^{or t^ *^o^o solution of 

the equation 


f^(l+t^)’^ + (l-02cos^9/t2)in2 , -[dln(f/Y)/dY](v./2vs1n^9) = — ^ . (12) 

sin'^e 

In eq. (11) 


2m 


- 2 . 


X 

where 


2dY^ 


Yod^ln(f/Y) 

dln(f/Y) 

2dY^ 

Y,Yq 

dlny 


W(t^.e),03) 


Y.Y, 


-3/2 

„ _ (3/2+sin^9/2)+sin^9]t’^(l+t^) - (cotVY2)ln2^^^ 


^ 0 ^ 0+tJ) ' + (l-s2cosS/t^)ln2 


-H 


2^ /i.2» 


( 14 ) 


max 
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3) Asympthotic Limits , These equations describe the synchrotron 
enlsslvlty for a general electron spectrum and pitch angle distribution. 

Before we can express j^(e) explicitly in terms of v and 9 we must 
specify the distribution function Below we shall consider three 

kinds of distribution. Before doing so we first describe some general 
features and comment on the accuracy of these expressions by considering some 
asympthotic limits. This will allow us to simplify these complicated expres- 
sions considerably. 

(a) High frequency, high energy limit. At high frequencies, at angles e 
such that sine Is not negligible, and for distributions so that 
\)^dln(f/Y) 

e « >0 « the right side of eq. (12) is much smaller than unity. 

2vdlny 

2 

This means that t^ » 1. In this limit the left hand side Is approximately 
equal to (l + 2^iLl1n_9.\ ^ |.q first order In e the solution 

3r \ lor / 

of eqs. (12) and (14) becomes 


t^/s1n^9 * e^yS 
0 ' 0 0 


2 

3eSln6 


» 1 . 


W 


3 


(15) 


Note that. In general, e may be a function of y (see examples given below), 
and that as 9 decreases y^ increases, which means that the largest contri- 

bution to the Integral comes from higher energy particles. Eq. (15) Is valid for 


1 a sine >> slnSj. » , (16) 

however, even for e = 9 (i.e. t^''» 1) it is accurate to better than 20 
percent. 
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(b) Sn»ll angle reglfue. For 0 < 9. ihe right hand side of eq. (12) 

c 

becomes larger than one so that t^ becomes less than unity, but Yq is still 

2 

much larger than one. In the extreme limit of 0 « 9^,, we have e/sin 9 » 1 

and t « 1. In this case the left hand side of eq. (12) is approximately 
1 2 

equal to "jf In so that for 0 « 0J, we have 
t 

t^/sin^e “ 02y2 In(sin9j./sin9), W* 1 . (17) 

Note that this equation will not be of much practical use since for e « 1 

eq. (15) is valid at all angles except for emission along the field line, 

where the emissivity is negligible (emissivity is zero at 0 ■ 0 except at 
the gyro- frequency). Thus, £v:cept for the unrealistic (i.e. not attainable 
in astrophysical situation) case of uniform magnetic field within the observed 
beam, we can ignore emission at angles 0 s 0_ and use eq. (15) (cf. Epstein 1973). 

W 

(c) Non-relativistic limit and low harmonics. At lower harmonics 
and/or for very steep electron spectra (which means primarily non-relativistic 
particles) — may exceed v/vj^ and the right side of eq. (12) becomes 
greater than one even at sine - 1. In the extreme case of e » 1, it is easy 

to see that, independent of e, the largest contribution to the emissivity 

comes from particles with non-relativistic energies. Thus, eqs, (14) and (12) 
reduce to 


g2y2 . 1/w » 1/e , e » 1. (18) 

(d) Test of the accuracy of approximations used in previous section. 

In the next section we will discuss the accuracy of the above approximate 
expression for three different kinds of distribution f(p,Y)» Before doing so 
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let us first examine the accuracy of the approximations used In the derivation 
of eqs. (11) to (14). In eq, (7) we assumed to be much less than unity. 
Clearly this wilt not be the case at all energies and angles. Our results 
would, however, be valid so long as this Inequality Is satisfied for the 
pitch angles and energies of particles which have the largest contribution 
to the Integral, namely, m » 3coso and y • y . For e « 1 and y 
given by eqs. (15) and (17), the expression / (1-z^) In eq. (7) 

Is equal to es1n/2 and c e^, respectively. Ignoring this term then the 
approximation » 0cos9 is justified^ as long as 

d1nf/dp < 2vcos9/vjjSln 8 . (19) 


v^dlnf/dp 

Note that at 9 ■ ,/Z eq. (7) Indicates that - 0 as long 

as dlnf/dy « v/vj^. As 9 decreases the Inequality 1n (19) becomes less 
restrictive. 


For e » 1 and Yq given by eq. (1C) Is no longer negligible, 

but In this case * gcosec^ * In ^ dlnf/dy » 0 as long 

as dlnf/dy << 2v/vjj. Thus, for particles with slowly varying pitch angle 

distribution, the relation * pcose is a good approximation throughout. 

In connection with eq. (4) It was mentioned that the approximation used 

2 3/2 

for the Bessel! function and their derivatives is valid if m(l-z ) * 

2 2 
■— — » 1- e » 1 from eq. (18) we have Yq * 1» « 1 

''b^o 

so that this condition is valid if v/vj^ s> 1, which is our basic assumption. 


n 


For e« 1 and sine 5 s slne^ from eq. (15), we find that the above 
Inequality reduces to 


(cslne)^'^ 


» 


Y*Y, 


1. 


( 20 ) 


which Is satisfied for most practical spectra. 

In the extreme limit of small angles e « e. the approximation for Bessel! 
function may not be valid since In this case the right side of the Inequality 
(19) Is replaced by In(s1ne./s1ne). However, this will happen at very small 

V 

angles which, as mentioned In connection with eq. (17), is of no consequence 
except at ultrarelativlstlc energies (see Epstein 1973, and references cited 
therein). 

Having justified all of our approximations, we can now substitute 
from eqs. (15) or (18) in eq. (11) to obtain to first order in e or 1/e 


p1n’®0exp{dln(f/Y)/dlnY^| ),e«l 


j^(9) - {.e\/c){vh^pf{$ cosB,y^n\ 2v/vJ° 

(ea slne/2) ^ e»l 
sln'^e ° 


. ( 21 ) 


As evident from expressions (13) and (14), the quantity X Is a compli- 
cated function of y^ and sine. However, as we have seen above, W(to,9) 
has a very simple form in the two limits. It turns out that, except for 
6«9 , eq. (14) can be simplified considerably. Combining the two asymthotic 
forms of W given In eqs. (18) and (15), v«! find that we can replace eq. (14) by 


W • 3/2 t 1 /(y| - 1) 


( 22 ) 


without loss of much accuracy. On Figure 1 we show the variation of the ratio 
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Of the actual value of W (from eq. 14) to Its approximate value given by 
eq. (22) with and sine. This shows that, except for the uninteresting 

case of sine « 1, the above expression Is an excellent approximation. 

We shall use eq. (22) Instead of eq. (14) In the following section. 
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II L SOME SPECIAL EXAMPLES 


In this section we derive the emissivity for three special particle 
distributions using the general eqs. (11 )» (12) and (21). Whenever possible 
we give a simple approximation to eq. (12). 

1) Thermal spectrum . We first consider the emissivity of thermal, non 
degenerate gas which was first investigated by Trubnikov (1961) and 
Drummond and Rosenbluth (1963) and which. has found application in fusion 


2 

There was some disagreement initially between Trubnikov and Drummond and 
Rosenbluth (1960) which was settled in favor of Trubnikov. Our results 
confirm this. 


plasma and in many astrophysical conditions such as solar flares, stellar 

coronae and accretion on white dwarfs (see e.g. Lamb and Masters 1979; Ramaty 

and Lingenfelter 1967; Ramaty and Petrosian 1972; Petrosian 1981). 

If the particle distribution is isotropic and Maxwellian with temperature 
2 

kT (in units of me ), then, in general 


- X.- 

f(lJ»Y)/Y “ Ce (y*-1)^ , 


where for kT < 1 


->5 1 


r -n I 

C « J 27r(kT)^ e^ (1 - 


3kT V 
16 


(23a) 


(23b) 


and n is the number density of the particles. From this we find that 



dy2 


+ 1 


<nn(f/Y)/<lY 


9 


(24) 
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Our detailed calculation with eq. (12) (described below) shows that the 
second term in the square brackets is always less than V|^/2v « 1 unless 
kT 2 1. However, if kT » 1, then the method of steepest descent used in the 
derivation of eq. (11) does not give a good approximation to the integral over 
particle ener:«ies. This is not of much concern because in most astrophyscial 
conditions (except in the early phase of the universe) such high temperatures 
are not encountered. Even when such temperatures are encountered, other high 
energy processes involving the numerous electron -positron pairs become more 
important than synchrotron radiation. Thus, we shall restrict ourselves to 
kT s 1, and in what follows immediately we ignore the square brackets in 
eq. (24). 

Let us first consider a Maxwellian gas with semi -relativistic temperature 
so that vkT/Ojj » 1. In this case e in eq. (12) is small (except for the 
uninteresting case of sine « Vj^/vkT) so that from eq. (15) we find 

3 ^* 1 , * 4vkT/3vjjSine, W = 3/2, X= (2kT/3Y^J^ (25) 

substitution of which in eq. (11) (keeping only the highest order term in y^)* 
or in eq. (21) gives 

1/3 

C(.kT/v^)exp{-^ [^(^)] }. (26) 

This is in agreement with results we obtain from Trubnikov's expression in this 
limit. In particular, for e = tt/2-(x with a « 1, the angular dependence 


reduces to 




„„ / v (>2 r 1 


} 


as given by Trubnikov (see also 


Bekefi 1966, p. 205). 
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For a non-relativistic gas with kT « 1 so that vkT/vjj is also much 
less than unity, the right hand side of eq. (12) is equal to Vjj/ 2 vkTsin 0 » 1 
so that its solution reduces to the simple expression in eq. (18): 

. 1 , . 2vI(T/y^. X • (2v/V|,)’‘ . (27) 

Substitution of this in eq. (11) or (21) gives 

j^(9) • (2 /?rre2vjj/c)C(vkT/vj,)^^^ ^ljS^jexp|-(v/vjj)ln(2vjj/evkTsin^9)|.(28) 

We find that the expression given by Trubnikov reduces to eq. (28) in this 
limiting case also. 

It is interesting that we find this agreement between the result here 
and those of Trubnikov even though the method of integration used here is 
quite different than that used by him and that eq. (12) bears no resemblance 
to his eq ’ation which plays the same role as our eq. (12). Presumably this 
agreement persists throughout. Note, however, that our result is more general 
and is applicable to a variety of spectral and pitch angle distributions 
(within the limitation specified in the previous section) while Trubnikov's 
result is applicable only to an isotropic Maxwellian gas with kT « 1. 

In the intermediate region, vkT/vj^ « 1, simple approximations such as those 
in eqs. (25) and (27) are not possible and we must solve eq. (12) numerically. In 
Figure 2 we show the variation of (y- 1) with vkT/vj^ for various values of 
sin 9 and for the limiting case of kT -► 0 (so that the square brackets in 
eq. (24) can be set equal to unity) and for kT = 1 (where we have included the 
square bracket in eq. (24) in our calculations). For kT -► 0 these curves 
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show a smooth transition between the two limiting expressions (25) and 
(27) except for very small value of s1n0. For kT « 1» vkT/vjj is always 
greater than unity and eq. (25} opears to be a good approximation throughout. 

From Figure 2 and eqs, (22) and (24) we find that the complicated eqs. (12) and be 
replaced by the following approximate equations without loss of much accuracy: 

[(2vkT/M|,)(l + 4.5vkTs1n^6/\»lj)''^’, kT « 1 
\ (4vkT/vi|,s1n8)^^^ , kT • 1 

(29) 

(2kT/Y*)(Yj- 1 )/(JyJ - 1) . kT s 1 . 


(yJ - 1) = 
x2 . 


These expressions asympthotically approach the expressions given in eqs. (25) 
ana (27). In addition, as shown by the inset in Figure 2 where we plot the 
ratio of - 1 fi^om (29) to that obtained from exact expression (12), 
the above approximation agrees with the exact result to better than few percent 
for v/vjj >> 1 and sine - 1 and to within a factor of 2 throughout. The 
largest deviation occurs at vkT/vj^ = 1 and sine « 1 , where eq. (17) is 
the appropriate asymthotic limit. Although we can improve on this agreement 
by adding new terms to eq. (29), the resultant improvement is not worth the 
sacrifice in the simplicity of eq. (29). 

To summarize, we can calculate the synchrotron emissivity from an isotropic 
Maxwellian gas at any kT ^ 1 using eqs. (11) and (29). Figures 3 
compares the results obtained from these equations with results from detailed 
integration of eq. (1). As is evident, these simple analytic expressions give 
excellent results not only at » 1, but surprisingly at lower harmonics 
where the accuracy of the steepest descent method is lower. 
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2) Power law energy spectrum . Power law spectra are other commonly 
used spectra In astrophysical problems. Usually power law spectra are defined 
with a low energy cutoff. To avoid such discontinuities and to simplify our 
analysis, we assume a spectrum of the form 

f(u.Y) “ 9 (m)c(i + (30) 

which converges at low energies; 

+1 

g(p)du - 1. (31) 

In eq. (29) e. plays the role of the low energy cutoff. For energies much 

greater than s the spectrum Is a power law with Index -« but tends to a con- 

stant at lower energies. The particles can be classified as ultra-relativistic 

or non-relativistlc If e » 1 or e « 1. We are interested primarily 

c c 

in cases with e_ » 1. For e. » 1 the maximum In the function 
c c 

used In derivation of eq. (11) becomes very broad so that the steepest descent 
method does not provide a good approximation to the Integrals. In this case 
one must use the well known ultra-relativistic results. 

For distributions which are not highly anisotropic (1.e.,dlng(p)/du « v/vj^ ) 
we can carry out calculation similar to that for a thermal gas. From eq. (30) 
we find 
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If vjjfi/v « 1, then » 1 so that eqs. (13), (14), (15) and (22) give 


®o * '• ’^0 " 3v(,(U«)s(ne • * * <’ > 


( 33 ) 


substitution of which in eqs. (11) or (21) gives 

/ f+i \ ( <s-l ) 

/ . , r- 1 -^^W3eX(«+1)sine\T' 

j^(^) • (/l7reSjj/2c) Cs^g(3^cos0)sine e ^ ^ 1 


(34) 


This expression has exactly the same dependence on \i, and e as the ultra- 
relativistic form (cf.e.g. Ginzberg and Syrovatskii 1964, p. 67). Even 
the numerical coefficients (which depend primarily on 6) agree to within a 
factor of 2, although the particle spectrum used in the ultra-relativistic 
calculations is different than that given by eq. (30). 

In the other extreme limit, V|^(1 +y)/\i >> 1, which is not of much interest 
here (see below) unless e. « 1, from eq. (18) we find that y « 1 and 

8* - -■p— . This is identical to eq. (27) if we identify as an equivalent kT. 

for the power law spectrum. In fact, substitution of this in eq. (11) or (21) 

gives an emissivity also identical to that in eq. (28) if we replace kT by 

-v/v. 6 

ej6 and set terms such as e °(1 + -r^) » 1. 

C 

In the intermediate ranges the expressions for y^ or 3^ are more com- 
plicated and depend on the value of e^.. Unlike in the case of thermal 
spectrum, it is not possible to give a simple expression for y^ - l because 
here y^ - 1 varies rapidly between its two asympthotic values: 
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y 2 . 1 . ^ Vk(1+6)/v « 1 (35a) 

° 3vjj(l+6)s1ne ° 

2vc_ 

2,1. £ V. (l+6)/\) » 1 (35b) 

in a manner similar to that found In degenerate Fermi distribution. 

In Figures 4 and 5 we plot (yJ - l)3vjj/(l+6)s1ne/4v vs v/5vjj 

for various values of e , sine and 6. We have attempted to fit the transi- 
tion region with an emperical formula. We have not found a simple expression 
which could reproduce the transition region even within a factor of 2 at all 
0 and v/v|j. However, as evident from Figure 4 for - 1, the high frequency 
form for yj - 1 (eq. 35a) provides a good approximation at all frequencies. 

On Figure 6 we compare the emissivlty obtained using eqs. (11), (12), (22), 

(32) and (35a) with results from numerical Integration of eqs. (1) and (2) 
for an Isotropic power law distribution with 6 = 3, * 1 and at e » 45°. 

Just as for the previous case, here again we find excellent agreement between spec- 
tra down to the gyro-frequency. Even the absolute value of the emissivlties at 
high V agree within 30 percent, This slight inaccuracy is because the lefthand 
side of expression (20) Is equal to 4, making this Inequality not a strong one. 

3) Non-isotropic pitch angle distribution . The two examples given above 
are valid if the inequality (19) is satisfied. This condition Is violated (specially) 
at low frequencies) if the particles have a highly anisotropic pitch angle distri- 
bution. For example, if f(p,Y) “ (vi-p^)° and if n > 2v/vjj, then the method of 
steepest descent cannot be used in the derivation leading to eq. (6). Neverthe- 
less, for a given n one can always find high enough frequencies where the above 
relations are valid. However, for a large n this could occur at such high 
harmonics where the emissivity is negligible. 
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On the other hand, as mentioned above, because of self-absorption the 
flux of radiation at lower harmonics Is negligible. Assuming that 
> 1 for V < 5v|j , then the Inequality (19) would be violated Ilf 
n a 10 for 9 » tt/2 and for even greater n at 6 < ir/2 (e.g., n > 20 at 

9 ■ n/4). It Is difficult to envision astrophysl cal conditions where the 
pitch angle distribution of the particles would be so highly anisotropic 
except If this anisotropy Is along the magnetic fldld lines, I.e. If particles 
form a tight beam along the field line. Even If particles are Initially 
accelerated along the field lines, their Interaction with the medium (through 
coulomb collisions, scattering by waves or by Inhomogeneities In the magnetic 
field) could quickly disperse the beam In the pitch angle space so that the 
Inequality (19) Is satisfied. However, if collisions are infrequent and the 
magnetic field Is very uniform, then the particle pitch angle could remain small. 

To complete our discussion let us consider a highly anisotropic distri- 
bution In pitch angle of the form 

f(|i,v) = (2/a2)e “o^“o f(y)* = 1 - ^ 


so that the inequality (19) is no longer satisfied. In this case the Taylor 
series expansion of eq. (5) gives 


j„(e) 


■?b" 


c 


(^)| dYfW) Yd.y.e) 

'‘»Jl 



2ada 




( 37 ) 


m 


= ^ (1-3COS6) and B = 1 + T-iMofe 

% '^b L4(i-0cose) 


+ 3coseln 


2(1-BC0S9) 

2 

3asin 9 - 
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As we shall see below, the largest contribution to the integral comes from elec 
trons with y of the order unity so that vya^/vjj « 1. Then we can neglect 
the slow logarithmic dependence of B in a and integrate eq. (37) to obtain 
to lowest order in vaJ/vj^ 




c 




es^sin^eyvo^ 

0 

4vjj(l ecosa) 


m 


f(y)dY . 


(38) 


Now following our previous procedure, we find that the largest contribu- 
tion to the integral comes from particles with 


» cose + 0 [l/ln(Vva|)] for 0 f ir/2 


3 


2 2. 


° ln(vjj/a2v) 


for 0 » ir/2 . 


Substitution of these in eq. (38) gives 


when 




4sin0v. 


A » cot exp In j ° 


'^b ecos^Oa^ 

0 { 


e v^n 

f(vj A. 
b 


V 1 / 'ecos^ea^ /-* 


A = *^ exp 


* '^b ea^v ^ I 


(4vjj/e^va2) , 


(39) 


(40) 


0 f w/2 , 


(41) 


0 = tr/2 . 


These are valid for semi -relativistic particles (i.e. for v/S\>^ or vkT/v^^ > 1). 
The situation becomes more complicated for non-relativistic particles since 
there exists two small parameters and vkT/v^j so that simple general 

expressions such as (40) are not possible. 
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IV. SUMMARY AND CONCLUSION 

We have found a single method of Integrating the complex synchrotron 
emissivlty over pitch angles and energies of ensemble of particles, 

1) Eqs, (11) to (14) present the results in their most general form. 

These equations can be used for evaluation of the synchrotron emissivlty 
from any particle distribution which falls fairly rapidly with energy and 
satisfies the inequality (19). The latter condition restricts the particle 
pitch angle distribution. 

2) In most practical applications (high frequencies and at directions 
away from the direction of the field lines), eqs. (12) and (21) provide a 
simple expression for the emissivlty. 

3) Our result when applied to an isotropic thermal gas of temperature T 
agrees, in two asympthotic limits vkT/v^j « 1 and » 1, with the earlier 
results of Trubnikov (1961) even though the method used here is quite different 
from Trubnikov* s, which is restricted tc emissivlty of a thermal gas. The 
simple emperical formula, eq. (29), connecting the two asympthotic limits is 
found to be sufficiently accurate for most practical application. 

4) Our result for a power law particle spectrum when extended to high fre 
quencies gives a result identical to that obtained from ultra-relativistic 
approximation as for its dependence on frequency, angle with respect to the 
magnetic field directions and on the parameters describing the particle 
distribution. 

5) For isotropic pitch angle distributions, spectra from our simple 
analytic expressions agree with numerical results to a high degree of accuracy 
even at lower harmonics, making resort to such detailed numerical calculations 
unnecessary. The absolute value of the synchrotron emissivity also agrees with 
20 to 30 percent at high frequency. 



23 


6) We also briefly describe a method for derivation of emlsslvlty from 
particles with highly anisotropic pitch angle distribution. Polarization and 
emission at lower harmonics and optically thick regime will be discussed In a 
subsequent paper. 
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FIGURii CAPTIONS 


Figure 1. Ratio of exact value of W (eq. IS) to Its approximate value 
(eq. 22) versus and o. Note that the ratio deviates from unity 
significantly when (Yq - 1) Is of order unity and only for sine « 1. 

For most practical purposes eq. (22) will be adequate. 

Figure 2. Variation of Yq - 1 (read Yq Instead of y In this figure) with 

vkT/vjj for an Isotropic Maxwellian gas at temperature kT. For sine of 

order unity Yq - 1 varies smoothly and monotonically between the two 
asympthotic values. Deviations from this behavior occurs at very small 
values of sine. Note that at high frequencies Is Independent of 
value of kT and depends only on vkT/vjj. Note kT ■ 0 means kT « 1. 

Inset. Ratio of the approximate value of yJ - 1 from eq. (29) to 

Its exact value given by eq. (12) for kT « 1 and kT h- 0. 

Figure 3. The synchrotron emissivlty normalized to the Raylelgh-Jean value 
((jiy * at e » tt/ 2 for kJ ^ 0.1 (or 50 keV for an electron). 

Points from our analytic equations using approximations In eqs. (11) and 
(29). Solid line fron numerical Integration of eq. (1) taken from Bekefi 
(1966) figure 6.10a, p. 203. Note the excellent agreement to much lower 
frequencies than expected from the analytic approximation. 

Figure 4. Deviation of yJ - 1 (y In the ordinate should be Yq) fnom Its high 
frequency asympthotic value versus v/6V|j at 0 * tt/2 for a power law 
spectrum. Note that for of order unity the asympthotic value provides 
an excellent approximation throughout. For « li Y«* 1 changes rapidly 

V w 

between the high and low frequency asympthotic values (eq. 35), with the 
transition value of v/6Vjj being Independent of 6 and e^. 
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Figure 5. Sanie as Figure 4 but for various values of sin9. Note the transition 

1 

becomes sharper with decreasing values of sin9 but occurs at the same value 
of v/Vjj6 for all 9. 

Figure 6* The synchrotron emissivity from an isotropic and power law particle 
distribution with 6 ■ 3 and 1 at 9 ■ 45°. Points from our 
analytic eq. (11) using the approximation in eqs. (22) and (35a). The 
solid line from the numerical integration of eq. (1) taken from Ramaty and 
Petrosian (1972). Note the excellent agreement to much lower frequencies 
than expected from the analytic approximations. 
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